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1. Introduction 

A many-body theory approach developed by the authors (Gribakin and Ludlow 2004) 
takes into account all main correlation effects in positron-atom interactions. These are 
polarization of the atomic system by the positron, virtual positronium formation and 
enhancement of the electron-positron contact density due to their Coulomb attraction. 
The first two effects are crucial for an accurate description of positron-atom scattering 
and the third one is very important for calculating positron annihilation rates. In this 
paper we apply our method to calculate the energies and annihilation rates for positron 
bound states with the hydrogen and halogen negative ions. 

Positron bound states can have a strong effect on positron annihilation with matter. 
For example, positron bound states with molecules give rise to vibrational Feshbach 
resonances, which leads to a strong enhancement of the positron annihilation rates in 
many polyatomic molecular gases (Gribakin 2000, 2001, Gilbert et al 2002, Gribakin et 
al 2010). Nevertheless, the question of positron binding with neutral atoms or molecules 
has been answered in the affirmative only recently, and the ubiquity of such states is 
only becoming clear now (Dzuba et al 1995, Ryzhikh and Mitroy 1997, Mitroy et al 
2001, 2002, Danielson 2009). On the other hand, it has been known for many decades 
that positrons bind to negative ions. Beyond the electron-positron bound state, or 
positronium (Ps), the simplest atomic system capable of binding the positron is the 
negative hydrogen ion (Ore 1951). The binding energy, annihilation rate and structure 
of the resulting compound, positronium hydride (PsH) have now been calculated to 
very high precision, e.g., by variational methods (Frolov et al 1997). The information 
available for heavier negative ions with many valence electrons is not nearly as accurate 
(see Schrader and Moxom 2001 for a useful review). Positronium halides have received 
most of the attention, with some calculations dating back fifty years; see, e.g., the 
Hartree-Fock calculations of Simons (1953) and Cade and Farazdel (1977), quantum 
Monte-Carlo work by Schrader et al (1992a, 1993), and more recent configuration 
interaction results (Saito 1995, 2005, Saito and Hidao 1998, Miura and Saito 2003). 

At first glance, the physics of positron binding to negative ions is much simpler than 
that of positron binding to neutrals. The driving force here is the Coulomb attraction 
between the particles. However, in contrast with atoms, the electron binding energy 
in a negative ion (i.e., the electron affinity, EA) is always smaller than the binding 
energy of Ps, |£i s | ~ 6.8 eV. This means that the lowest dissociation threshold in the 
positron-anion system is that of the neutral atom and Ps, so that the positron bound 
to a negative ion may still escape by Ps emission. Hence, for the system to be truly 
bound, its energy should lie below the Ps-atom threshold, and the positron energy in 
the bound state, e , must satisfy 

|e | > \E ls \ - EA. (1) 

This situation is similar to positron binding to atoms with ionization potentials smaller 
than 6.8 eV. The structure of these bound states is characterized by a large "Ps cluster" 
component (Mitroy et al 2002). 
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The proximity of the Ps threshold in anions means that to yield accurate binding 
energies, the method used must account for virtual Ps formation. The positron 
also polarizes the electron cloud, inducing an attractive polarization potential. It 
has the form — ae 2 /2r 4 at large positron-target separations r, where a is the dipole 
polarizability of the target. The method should thus be capable of describing many- 
electron correlation effects. The halogen negative ions have np 6 valence configurations, 
and are similar to noble-gas atoms. Here many-body theory may have an advantage 
over few-body methods (Dzuba et al 1996). 

The many-body theory approach developed by the authors (Gribakin and Ludlow 
2004), employs B-spline basis sets. This enables the sum of the electron-positron ladder 
diagram sequence, or vertex function, to be found exactly. This vertex function accounts 
for virtual Ps formation. It is incorporated into the diagrams for the positron correlation 
potential and correlation corrections to the electron-positron annihilation vertex. To 
ensure convergence with respect to the maximum orbital angular momentum of the 
intermediate electron and positron states in the diagrams, Z max , we use extrapolation. 
It is based on the known asymptotic behaviour of the energies and annihilation rates as 
functions of Z max (Gribakin and Ludlow 2002). 

In Gribakin and Ludlow (2004) the theory has been successfully applied to positron 
scattering and annihilation on hydrogen below the Ps formation threshold, where the 
present formalism is exact. This theory is now extended to treat the more difficult 
problem of positron interaction with multielectron atomic negative ions. We test the 
method for H~, and consider the halogen negative ions F~, Cl _ , Br - and I - . Note that 
a conventional notation for the positron bound with a negative ion A" is PsA, rather 
than e + A~, hence one has PsH, PsF, etc. (Schrader 1998). 

2. Calculation of positron binding using Dyson's equation 

The many-body theory method for positron bound states is similar to that developed 
for electron-atom binding in negative ions (Chernysheva et al 1988, Johnson et al 1989, 
Gribakin et al 1990, Dzuba et al 1994), and used for positron-atom bound states by 
Dzuba et al (1995). 

The Green function of the positron interacting with a many-electron system 
("target") satisfies the Dyson equation (see, e.g., Migdal 1967), 

(E - H )G E (r, r') - J S B (r, r")G E (r", r')dr" = 8(r - r'), (2) 

where H is the zeroth-order positron Hamiltonian, and T, E is the self-energy. A 
convenient choice of Hq is that of the positron moving in the field of the Hartree-Fock 
(HF) target ground state. The self-energy then describes the correlation interaction 
between the positron and the target beyond the static-field HF approximation (Bell 
and Squires 1959). It can be calculated by means of the many-body diagrammatic 
expansion in powers of the electron-positron and electron-electron Coulomb interactions 
(see below). 
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If the positron is capable of binding to the system, i.e., the target has a positive 
positron affinity PA, the positron Green function Ge(f, r') has a pole at E = Eq = —PA, 

GeM = (3) 

Here ipo( r ) is the quasiparticle wavefunction which describes the bound-state positron. 
It is equal to the projection of the total ground-state wavefunction of the positron and TV 
electrons, \l/o( r i> • • • > r N, r ), onto the target ground-state wavefunction, $o( r i, • • • > t n) , 



Mr) 



J $*(ri, . . . , rjv)* (ri, ■ . . , rjv, r)dr x . . . dr^. (4) 
The normalization integral for ipo(r), 

l^o(r)| 2 dr<l, (5) 



can be interpreted as the probability that the electronic subsystem of the positron-target 
complex remains in its ground state. 

The magnitude of a quantifies the extent to which the structure of the complex 
is that of the positron bound to the anion, e + A~, as opposed to that of the Ps atom 
orbiting the neutral atom. Such separation is a feature of the "heuristic wavefunction 
model" (Mitroy et al 2002). If the former component dominates the wavefunction, the 
value of a is expected to be be close to unity. If, on the other hand, the wavefunction 
contains a large "Ps cluster" component, the value of a will be notably smaller. 

By taking the limit E — » 6q in equation (j2J), one obtains the Dyson equation for 
the quasiparticle wavefunction ipo(r) ("Dyson orbital") and the bound-state energy 6q, 

ifoV'o(r) + J S eo (r, r')^ (r')dr' = e Q %{r). (6) 

This equation is analogous to the standard Schrodinger eigenvalue problem, except that 
the correlation potential E depends on the energy. 
The eigenstates <£> e (r) of the Hamiltonian H , 

H if £ {r) = e(p 6 (r), (7) 

are characterized by their energies e and the orbital angular momentum /, implicit in 
this notation. They form a complete single-particle positron basis set. For negative ions 
the spectrum of H consists of both discrete (Rydberg) and continuum states. This set 
can be used to expand the quasiparticle bound-state wavefunction, 

Vo(r)=^C £ ^(r), (8) 

e 

and cast the Dyson equation ([6]) in the matrix form, 

eC £ + Y J {^e \e')C £l =e C £ , (9) 

£' 

where the sums in equations (jEJ) and ([9]) include the positive-energy continuum, as well 
as the discrete negative-energy states, and 



{e\^ E \e') = J ^(r)£ £ (r,rOMOdrdr'. 



(10) 
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In practice the continuum is discretized by using a B-spline basis set (see below), 
and equation (Q becomes a matrix eigenvalue problem. Therefore, to find Eq and C e 
one simply needs to diagonalize the matrix, 

s £ 5 e£l + (e\E E \s'). (11) 

Its lowest eigenvalue s (E) depends on the energy E at which T, E is calculated, and 
the diagonalization must be repeated several times until self-consistency is achieved: 
Eq(E) = E. Knowing the dependence of the eigenvalue on E allows one to determine 
the normalization integral, equation (j5]), via the relation (Migdal 1967), 

de (E)' 



dE 



(12) 



E=e , 



Note that owing to the spherical symmetry of the target, the states and ipo have 
definite orbital angular momenta. To find the bound positron ground state it is sufficient 
to calculate the self-energy matrix ([TO]) for the s-wave positron only. 

The accuracy of the binding energy obtained from the Dyson equation depends 
upon the accuracy to which the self-energy has been determined. As mentioned in 
the Introduction, polarization of the target and virtual Ps formation are the two most 
important effects that need to be accounted for. The effect of target polarization is 
described in the leading order by the 2nd-order diagram E< 2 \ fi gure [l](a). Following 
Gribakin and Ludlow (2004), the Ps formation contribution £( r ) shown in figure [T^b), is 
obtained by summing the electron-positron ladder diagram series to all orders, figure [2j 
This procedure amounts to calculation of the electron-positron vertex function T shown 
by the shaded block. Analytical expressions for the diagrams can be found in Gribakin 
and Ludlow (2004). 




n 

(a) 




Figure 1. Main contributions to the positron self-energy. Diagram (a) describes the 
effect of polarization in the lowest, second order, £( 2 ); diagram (b) accounts for virtual 
Ps formation, Y,( r \ Top lines in the diagrams describe the positron. Other lines with 
the arrows to the right are excited electron states, and to the left - holes, i.e., electron 
states occupied in the target ground state. Wavy lines are the Coulomb interactions. 
Summation over the intermediate positron, electron and hole states is carried out. 

For hydrogen, the self-energy is given exactly by the two diagrams in figure [U 
£ = £( 2 ) + S^ r \ provided the intermediate electron and positron states are calculated in 
the field of the bare nucleus (Gribakin and Ludlow 2004). For multi-electron targets one 
may also consider higher-order corrections not included in the virtual Ps contribution, 
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Figure 2. Electron-positron ladder diagram series and its sum, the vertex function T 
(shaded block) . Comparison between the left- and right-hand sides of the diagrammatic 
equation shows that T can be found by solving a linear matrix equation. 

figure Wfo)- In the present calculation a set of 3rd-order diagrams shown in figure [3] 
will be included. Diagrams (a), (b), (c) and (d) represent corrections to the 2nd-order 
polarization diagram, of the type described by the random-phase approximation with 
exchange (RPAE, Amusia et al 1975). They account for the electron- hole interaction and 
screening of the positron Coulomb field, and correct the value of the dipole polarizability 
a of the target. Diagram (e) accounts for the positron-hole repulsion. The contribution 
of the diagrams in figure |3] is denoted collectively by . Calculation of these diagrams 
will allow us to gauge the importance of these corrections, and even to include effectively 
higher-order diagrams (see below). 




(a) (b) (c) (d) (e) 



Figure 3. Third-order contributions to the positron self-energy, E^ 3 ^. Mirror images 
of the diagrams (c) and (d) are also included. The top line describes the positron. 

Of course, any many-body theory calculation can at best include only dominant 
classes of diagrams, leaving out an infinite number of other higher-order diagrams. 
For example, the diagram in figure H] has the effect of screening the positron-electron 
interaction accompanying virtual Ps formation in TS T \ figure [T^b). 



£' Vj V 2 £ 




Figure 4. Screening correction to the virtual Ps contribution E^ r ). 

The number of diagrams increases rapidly as one moves to higher orders. The 
effort required to evaluate these diagrams becomes prohibitive relative to their small 
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contribution. It would therefore be useful to find a simple method to estimate the 
contribution of the higher-order diagrams. This will allow us to take into account the 
effect of electron screening beyond the corrections shown in figure [3j 

A useful quantity for estimating the size of a contribution to the self-energy E is a 
dimensionless measure of its "strength" (Dzuba et al 1994), 

<te(£)= fG(°\r', r)Z E (r,r')drdr' (13) 
= _J2 {£ ^ eK) , (14) 

V 

where is the Oth order positron Green function calculated at E = 0. Let 

Se - g-jm (15) 

be the ratio of the strength of the sum of the 3rd-order polarization diagrams (figure [3]) 
to the strength of the 2nd-order polarization diagram. This quantity Se can then be 
used to estimate higher-order contributions to £( r ). 

As a check, we test that the binding energy obtained using £( 2+3 ) = E^ 2 ) + S^ 3 ^ is 
close to that obtained with T,^ multiplied by 1 + Se, i.e., using 

(l + ^)£( 2 ) =£( 2+3 '). (16) 

An estimate of the total self-energy corrected for the screening effects in the lowest order 
can then be obtained as 

(1 + S E ) [X (2) + £ (r) ] • (17) 

As we will see in section HJ the relative effect of screening is negative, Se < 0, 
which means that screening reduces the magnitude of the self-energy. Similar higher- 
order terms in the self-energy expansion will alternate in sign. Thus, for example, 
a diagram such as that shown in figure |5l will tend to compensate the lowest-order 
screening correction in figure HI 

£' £ 




Figure 5. Higher-order screening correction to the virtual Ps contribution £( r ), cf. 
figure HI 

Assuming that the sequence of screening corrections behaves like a geometric series, 
its effect can approximately be taken into account by using the screened self-energy, 
which we denote S^ scr ^ 

S (scr) = _J_ r E (2) + S (D1 (lg) 

1 — Se 



Positron binding to negative ions 



8 



We expect that this approximation should yield our best prediction for the binding 
energy. 



3. Calculation of the positron annihilation rate in the bound state 

The spin-averaged positron annihilation rate T a in the bound state can be expressed in 
terms of the average contact electron-positron density p ep (see, e.g., Berestetskii et al 
1982) 

T a = nr^cpep, (19) 

where r is the classical radius of the electron, c is the speed of light, and p ep is given 
by the integral, 

N f 

Pe P = Yl / l^o(ri,r 2 ,...,r iV ,r)| 2 5(r-r i )dr 1 ...dr Ar dr, (20) 

8=1 

where \l/o( r i; r 2, • • • , r N, r ) is the full (N + l)-particle bound-state wavefunction of the 
TV electron coordinates and positron coordinate r. 

Figure [6] shows a series of diagrams that would constitute a complete set of 
annihilation diagrams for a one-electron system (Gribakin and Ludlow 2004). Here e 
represents the positron bound-state Dyson orbital ^>o( r ) normalized as per equation (jSJ). 




(e) 



(f) 



Figure 6. Many-body theory expansion for the contact electron-positron density. The 
solid circle in the diagrams is the delta- function annihilation vertex, cf . equation (|20[) . 
Diagrams (b), (d) and (e) are multiplied by two to account for their mirror images. 

Diagram (a) in figure [6] is the overlap of the positron and HF electron densities. It 
represents the independent-particle approximation to the annihilation vertex with the 
contact density 



N 



JO) 
rep 



n=l 



£ / |0 n (r)|Vo(r)| 2 dr, 



(21) 



where n (r) is the HF orbital of hole n. The first-order correction, diagram (b), can be 
thought of as the analogue of £( 2 ), and will be denoted by pip . The diagrams with the 
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vertex function T, e.g., (d) in figure El are particularly important in the calculation of 
p ep , as the annihilation takes place at a point, and is strongly enhanced by the Coulomb 
attraction in the annihilating electron-positron pair. 

The diagrams shown in figure [6] represent a basic set of contributions with a single 
hole line, which one needs to consider to obtain a reliable answer. We will denote 
the corresponding density as pfp + pip \ pip^ representing a correction to the zeroth- 
order contact density. Similar to the self-energy diagrams in figure (H they represent 
the exact answer for the positron-hydrogen system, provided the electron and positron 
intermediate states are calculated in the field of the bare nucleus (Gribakin and Ludlow 
2004). For complex many-electron systems it may be necessary to account for the effects 
of electron screening when calculating p ep . A series of RPA-type annihilation diagrams 
is therefore also calculated, see figure UJ 




(a) (b) (c) (d) 




(e) (f) (g) 

Figure 7. Annihilation diagrams with two Coulomb interactions, including those of 

(2) 

RPA-type, p y ep ' . The top line describes the positron. All the diagrams have equal 
mirror images. 

Similarly, screening corrections to the annihilation diagrams containing the T block 
can also be considered, as shown in figure 




Figure 8. Screening correction to the annihilation diagram containing the vertex 
function T. 

The diagrams shown in figure [71 pip , can be thought of as next-order corrections 
to p e p , diagram (b) in figure El By evaluating the ratio 

C = P ( S/P < $, (22) 
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an estimate can be made of the total contact density p ep that includes higher-order 
corrections in a manner similar to the self-energy [cf. equation ffl8|) ]. 

1 

~c 



Pe P = Pe¥ + ^^P { e?- (23) 



4. Numerical Implementation 

The Hartree-Fock ground state of the negative ions is first found. The frozen- 
core HF Hamiltonian for an electron or a positron (with and without exchange, 
respectively) is then diagonalized in a B-spline basis (Sapirstein et al 1996). The 
corresponding eigenvectors provide bases of single-particle electron and positron states, 
cf. equation (jTJ). The spectrum of these states for the electron includes the negative- 
energy ground-state orbitals (hole states) and positive-energy excited states spanning 
the electron continuum, see, e.g., figure 6 in Gribakin and Ludlow (2004). The positron 
basis contains a number of negative-energy Rydberg states augmented by the discretized 
positive energy positron "continuum" . 

The effective spanning of the continuum is achieved by using an exponential radial 
knot sequence for the B-splines. For H~, the first 23 eigenstates generated from a 
set of 60 splines of order 9 were used, with a box size of R = 30 au. For the other 
systems, namely F~, CI - , Br~ and I - , the first 20 states from a set of 40 splines of 
order 6 were used with R = 30 au. Only the outermost s and p subshells were included 
when calculating the self-energy and annihilation diagrams. More strongly bound inner- 
shell electrons are only weakly perturbed by the positron. Their contribution to the 
correlation potential and annihilation vertex is relatively small, and has been neglected. 

The diagrammatic contributions to the self-energy and contact density described in 
sections [2] and [31 are calculated by direct summation over the intermediate electron and 
positron states, and the vertex function is found by solving a linear matrix equation. 
The use of B-spline bases ensures quick convergence with respect to the number of 
states with a particular angular momentum I included in the calculation. In addition, 
the convergence with respect to the maximum orbital angular momentum included in 
the calculation, Z max , also needs to be considered. This is done by extrapolation through 
the use of the asymptotic formulae (Gribakin and Ludlow 2002), 

£ - £lU -CTW (24) 



and 



B 



p " = »^' + (i^TWy (25) 

where 4' max) and Pep ^ are the bound-state energy and contact density obtained in 
a calculation for a given Z max , and A and B are constants. While the derivation 
of equations ( 1241 and ( 1251) is based on perturbation theory (Gribakin and Ludlow 
2002), this asymptotic behaviour is confirmed by nonperturbative many-body-theory 
calculations (Gribakin and Ludlow 2004), and by configuration-interaction calculations 
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of positron binding and annihilation in atoms (Mitroy and Bromley 2006). The constants 
A and B are found numerically and have different values for each system studied. 

5. Results 

5.1. Details of calculations for PsCl 

In this section a detailed examination of the results for PsCl will be presented. This 
should illustrate how the final results for the other systems were arrived at. 

The positron radial wavefunctions for PsCl obtained by solving the Dyson equation 
using the ah initio self-energy £( 2+r+3 ) and the screened self-energy £ < - scr - ) are compared 
to the HF positron wavefunction in figure |9j The inclusion of the attractive correlation 
potential results in lower energies of the positron bound states, and hence, more compact 
wavefunctions. Otherwise the two Dyson orbitals are very similar to each other. 

0.6 

o 

'to 0.4 

3 

Hi 
> 
CO 

3= 

co 

=5 0.2 
co 



"0 5 10 15 20 

Radius (au) 

Figure 9. Radial positron wavefunction in PsCl: , Hartree-Fock; — • — , Dyson 

orbital calculated using the self-energy £( 2 + r + 3 ); - — ; Dyson orbital obtained with 
^(scr)^ rpke HF and Dyson orbitals on the graph are normalized to unity. 

When solving the Dyson equation, the self-energy and the bound-state energy e 
were calculated for a number of maximum orbital angular momenta, e.g., / max = 7-10, 
and then the asymptotic behaviour (I24I) was used to find the result for Z max — » oo. This 
procedure is illustrated for PsCl in figure [TQJ Extrapolation from Z max = 10 to infinity 
increases the binding energy by about 0.5%. 

Before the contact density p ep can be determined, the positron Dyson orbital must 
be correctly normalized via equation ([121) . This is achieved by calculating the self-energy 
for a number of energies E and finding the lowest eigenvalue of the matrix ( fPTl) at these 
energies, giving £q(E). This is repeated to self-consistency, E = £q(E), and the gradient 
deo(E)/dE is found at this point. This is illustrated for PsCl in figure [TT1 
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Figure 10. Convergence of the binding energy for PsCl as a function of i max - Open 
circles connected by a solid line to guide the eye, show the energies £o calculated using 
^(2+r+3) a j. £ _ „o.207 au for Z max = 7-10; dashed line shows extrapolation. 




Figure 11. Open circles connected by the solid line show the positron energy £q(E) 
for PsCl obtained from the Dyson equation with I]( 2 + r + 3 ) ) as a function of the energy 
at which the self-energy was calculated; dashed line is £o = E. The intersection of the 
two lines, £o(E) — E, gives the binding energy. The gradient of £o{E) at this point is 
used to calculate the normalization constant a from equation ([5]). 



According to equation (125]) . the electron-positron contact density p ep converges 
much slower than the energy, and extrapolation with respect to Z max is much more 
important here. This is illustrated for PsCl in figure IT2l Extrapolation beyond Z max = 10 
increases the contact density, and hence, the annihilation rate, by about 30%. 

The calculations are performed for a number of approximations to the self-energy. 
This enables us to determine the relative magnitude of various diagrams, and helps 
to clarify which physical effects are important to include so as to obtain an accurate 
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Figure 12. Convergence of the electron-positron contact density p ep for PsCl as a 
function of Z max . Open circles connected by a solid line to guide the eye, show values 
of pfp + pip^ + pep obtained for Z max = 7-10 in the calculation using £( 2+r + 3 ); dashed 
line shows extrapolation. 

binding energy. The positron energies obtained for PsCl using different approximations 
are given in table [TJ 



Table 1. Positron bound-state energies £o (in a u) for CI obtained using various 
approximations to the correlation potential. 



HF 


E(2) 


E (2+r) 


£(2+3) 


s (2+r+3) 


£(2+3') £(scr) 


-0.1419 


-0.1855 


-0.2276 


-0.1663 


-0.2072 


-0.1641 -0.1998 



In the HF approximation, the energy of the lowest s wave positron state is 
—0.1419 au. The Ps binding energy (BE) is determined from the positron affinity 
PA = |eo| using equation ( )26|) . which gives BE = 0.672 eV. This value is in agreement 
with the HF results of Cade et al (1977). Therefore, PsCl is bound even in the static HF 
approximation. However, the self-energy is essential in determining an accurate binding 
energy. The second-order polarization diagram 

S (2) 

increases the binding energy, and 
the inclusion of the virtual Ps formation contribution £( r ) increases it even further. 
When we add the 3rd-order corrections to E^, the binding reduces noticeably. 
This means that screening of the Coulomb interaction is important. In particular, the 
positron binding energy of 0.2276 au, obtained with £( 2+r ), becomes equal to 0.2072 au 
when the total self-energy £( 2+r+3 ) is used. 

However, this calculation neglects the effect of screening on the virtual Ps-formation 
contribution Evaluating the magnitude of screening via equation ( TT5|) we obtain 

Se = —0.45. When the effect of screening on is included via the factor 1 + Se, 
as per equation ( fl6l) . the corresponding result (labelled £( 2+3 ')) is very close to that 
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obtained with E( 2+3 ) (see table [T]). The application of approximation (ITS]) , denoted 
gives our best estimate of the positron binding energy, 0.1998 au, corresponding 
to a Ps binding energy of 2.245 eV. This value is only slightly below the completely ab 
initio value of 2.437 eV obtained using £( 2+r + 3 ). 

The contact densities calculated using the Dyson orbitals obtained with £( 2 + r + 3 ) 
and £ <scr ) are quite close, as are the energies and wavef unctions. We show the 
breakdown of the contributions to p ep in table [2j To appreciate the scale of densities 
involved, it is useful to remember that the contact density of ground state Ps is 
p e p(Ps) = l/8?r « 0.0398 au. Note that although £( 2 + r + 3 ) gives a slightly larger binding 
energy and a more compact positron wavefunction than E^ scr \ the densities obtained in 
the former approximation are lower. This is due to a smaller normalization constant a, 
which results from a somewhat stronger energy dependence of £( 2+r + 3 ). 

Table 2. Breakdown of contributions to the electron-positron contact density in PsCl 
(in au). 



Approx. 


(0) 
Pep 


Pep 


(0) , (A) 
Pep ~r Pep 


p^ep total 


Pep, eq. (J53J) 


E (2+r+3) 


0.00841 


0.00931 


0.04263 


-0.00155 0.04108 




£(scr) 


0.00872 


0.00964 


0.04444 


-0.00162 0.04281 


0.03929 



The zeroth-order diagram, p ep , gives only about 20% of the total density, with p ep 
giving another 20% and the rest coming from higher order diagrams in pi p ^ (figure [6]). 
As with the self-energy, the inclusion of screening effects (pip , figure [7]) reduces the 
total. However, the effect of screening on the annihilation vertex is much smaller than 
that on the correlation potential, as indicated by the value of C — —0.17, equation f l22|) . 
Physically, this is related to the fact that in the annihilation vertex corrections, small 
electron-positron separations dominate. Finally, using equation (|23|) to account for the 
effect of screening on the diagrams in pip , we obtain our best prediction for the contact 
density (last column in table |2]). This corresponds to the PsCl decay rate of 1.984 ns -1 , 
which is close to the spin-averaged decay rate of Ps, 2.01 ns -1 . 

5.2. Results for PsH, PsF, PsCl, PsBr and PsI 

The final results for PsH, PsF, PsCl, PsBr and PsI obtained with the correlation 
potential £( 2 + r + 3 ) and density pip + pi^ + pf p \ and £( scr ) and screened densities from 
equation f[2"3"j) . are shown in table El In all cases the positron is bound in the s-wave, 
all higher lying quasi-bound states being unstable against Ps emission. Note that the 
latter is true for the electron-spin-singlet states, as excited "unnatural parity" electron- 
spin-triplet Ps-atom bound states have been discovered recently for the hydrogen and 
the alkalis (Mitroy and Bromley 2007, Mitroy et al 2007). 

The positron binding energy is highest in PsH. This is a consequence of the small size 
of the hydrogen atom, and the small value of its electron affinity, which makes for strong 
electron-positron correlation effects. Beyond PsH the binding energy decreases along 
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Table 3. Positron binding energies, normalization constants and contact densities for 
PsH and Positronium halides. 



Compound 


£o a au 




Pep 


£o b au 




Pep 


PsH 


-0.27619 


0.714 


0.05231 


-0.26338 


0.748 


0.05037 


PsF 


-0.22778 


0.950 


0.04790 


-0.22489 


0.958 


0.04913 


PsCl 


-0.20718 


0.875 


0.04108 


-0.19975 


0.894 


0.03929 


PsBr 


-0.20373 


0.834 


0.03910 


-0.19523 


0.868 


0.03788 


PsI 


-0.19805 


0.794 


0.03707 


-0.18878 


0.835 


0.03582 



a Dyson equation solved using £( 2 + r + 3 ) ; density p ep = pip + pip + Pep ■ 
b Dyson equation solved using J]( scr ) j density p ep from equation ([23| . 



the halogen sequence, mostly due to a stronger positron repulsion from the positively- 
charged atomic cores in heavier systems. 

Values of the normalization parameter a in table E] give some insight into the 
structure of these compounds. PsH has the smallest value of a and its structure therefore 
has a large component that describes Ps bound to the neutral atom ("Ps cluster"), the 
small electron affinity of H playing a role in this. PsF has the largest value of a and 
its structure can best be described as a positron bound to F~. Generally, all of the 
compounds considered have large values of a. This indicates that a positron bound to 
the negative ion is the dominant component of the structure. This is a consequence of the 
stable noble-gas-like structure of the halogen negative ions. In contrast, positron bound 
states with the weakly-bound alkali negative ions have a distinct Ps-atom character 
(Mitroy et al 2002). 

It is interesting to compare the positron wavefunctions obtained from the Dyson 
equation for H~ and the halogen anions. In figure [13] the wavefunctions obtained 
with the self-energy £( 2+r + 3 ) are shown. The shape of the positron wavefunction is 
determined by a balance between the Coulomb and correlation-potential attraction at 
large separations, and the Coulomb repulsion from the nucleus at smaller radii. The 
positron wavefunctions for PsH and PsF are quite similar. This feature reflects the high 
positron binding energy to the H~ and F~ ions, and the fact that the corresponding 
atoms have the smallest radii. As the positron binding energy decreases, the positron 
wavefunction relaxes outwards. This feature is seen as one moves along the halogen 
sequence. One can also observe the increasing "expulsion" of the positron from the 
atomic core region, caused by the Coulomb repulsion from the nucleus. 

5.3. Comparison with other theoretical results 

The positron binding energies can be converted to Ps binding energies via the simple 
relation, 

BE(PsA) = EA(A) + PA(A-) + E ls (Ps), (26) 

where BE(PsA) is the binding energy of a Ps atom to a generic atom denoted by 
A, EA(A) is the electron affinity of the atom, PA(v4~) is the positron affinity of the 




Figure 13. Comparison of the radial Dyson orbitals obtained with £( 2+r + 3 ): solid, 
PsH; dashed, PsF; long-dashed, PsCl; dot-dashed, PsBr; dot-double-dashed, PsI. For 
the purpose of comparison, all orbitals are normalized to unity. 



negative ion and £a s (Ps) = —6.8028 eV is the spin-averaged energy of Ps. For the 
electron affinities needed, the values of 0.7542 eV for H, 3.4012 eV for F, 3.6127 eV for 
CI, 3.3636 eV for Br and 3.0590 eV for I, have been used (Andersen et al 1999). 

Annihilation rates r a in units of ns _1 are obtained by dividing the contact density 
Pep by the conversion factor 10 9 s/(7rrQc) = 0.0198 ns x au. The final Ps binding 
energies and positron annihilation rates are shown in table HI and compared with other 
calculations and experiment. 

For PsH very accurate variational calculations are available (Frolov et al 1997). 
Our many-body theory calculations are in good agreement with these results, and both 
the Ps binding energy and positron annihilation rate are accurate to within 5%. Because 
of the small electron affinity of hydrogen, PsH is more difficult for many-body theory 
to treat than larger, more tightly bound systems with many valence electrons. The 
results for the heavier systems should therefore be of similar or possibly even of greater 
accuracy than the results for PsH. 

For PsF a few theoretical calculations are available, the present Ps binding energy 
and positron annihilation rate agreeing most closely with multi-reference configuration- 
interaction calculations (Miura and Saito 2003, Saito 2005). The Ps binding energy is 
also close to a diffusion Monte-Carlo calculation by Bressanini et al (1998). 

The present Ps binding energies for PsCl, PsBr and PsI are greater than those 
obtained using a second-order variational perturbation method (Saito et al 1995, 1998) 
and Monte-Carlo calculations (Schrader et al 1992a, 1993). Our values are in better 
agreement with, although consistently smaller than, multi-reference configuration- 
interaction calculations by Saito (2005). The positron annihilation rates from the present 
calculation and that of Saito (2005) are in reasonable agreement, although our values 
are consistently higher. 
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Table 4. Ps binding energies and positron annihilation rates for PsH and positronium 
halides compared with other calculations and experiment. 



Compound 


Present 


results 


Other 


results 




Ps BE (eV) 


r Q (ns- 1 ) 


Ps BE (eV) 


r a (ns- 1 ) 


PsH 


1.118 


2.544 


1.066126 a 
l.liO.2 1 


2.4361 a 




O 71 Q 

1. 1 lo 


z.4oz 


z.sUo 
2.70 c 
2.24 d 

i.yo±u.i / s 

2.838 f 

2.9±0.5 j 


z.Uly 
1.98 c 


PsCl 


2.245 


1.984 


2.350 b 
1.91± 0.16° 
1.62 d 
2.0±0.5 j 


1.504 b 


PsBr 


1.873 


1.913 


2.061 b 

1.14±0.11 g 

1.25 h 


1.371 b 


PsI 


1.393 


1.809 


1.714 b 
0.56 h 


1.254 b 



Theory: a Frolov et al 1997, b Saito 2005, c Miura and Saito 2003, d Saito 1995, 
c Schrader et al 1992a, f Bressanini et al 1998, s Schrader et al 1993, Saito et al 
1998. 

Experiment: 1 Schrader et al 1992b, j Tao et al 1969. 



5.4- Comparison with experiment 

For PsH, a direct experimental measurement of the Ps binding energy (Schrader et al 
1992b) is available. In this experiment, the reaction e + + CH4 — > CH^ + PsH was studied 
by detecting the CHj~ ions. From the experimentally determined threshold energy for 
CH3 production and the various bond energies, a value of 1.1 ± 0.2 eV was found for 
the PsH binding energy. This value is in excellent agreement with theory, though much 
less accurate. 

So far, there have been no direct experimental measurements of the Ps binding 
energy for the halogens, however estimates of the Ps binding energy for PsF and PsCl 
has been made (Tao et al 1969), see table HI The PsCl binding energy was estimated 
by studying positron annihilation in CI2 and Ar-Ck* gas mixtures. The appearance of a 
shoulder in the positron annihilation lifetime spectrum was attributed to the reaction, 
Ps+Ck — > PsCl + CI. From a knowledge of the energy at which this shoulder begins and 
the CI2 dissociation energy, a binding energy of about 2.0 eV was estimated for PsCl. 
The estimate of the PsF binding energy was obtained from the observation that when a 
hydrogen atom in benzene is replaced by fluorine, the fraction of positrons annihilating 
with the longest lifetime, as ortho-positronium, was reduced from 40 to 27%. This 
reduction was assumed to be due to the reaction, C 6 H 5 F + Ps — > PsF + C 6 H 5 . From 
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a knowledge of the threshold energy and the relevant dissociation energies, a binding 
energy of about 2.9 eV was estimated for PsF. The present results support these early 
estimates. 

6. Concluding remarks 

Traditionally, many-body theory has had more success in treating purely electronic 
systems (Chernysheva et al 1988, Dzuba and Gribakin 1994), than systems that contain 
a positron. In particular, earlier many-body theory calculations relied on very simple 
(Amusia et al 1976) or approximate (Dzuba et al 1995, 1996) treatments of the virtual 
Ps formation contribution to the correlation potential. However it is now clear that 
many-body theory is capable of giving positron and Ps binding energies and positron 
annihilation rates that are accurate to within a few percent for many-electron systems. 

The calculated binding energies and annihilation rates for positronium halides 
should serve as a useful reference for other theoretical calculations and future 
experiments. We have also performed an extended analysis of various contributions 
to the correlation potential and contact density, especially of the role of screening. This 
will be helpful for the problem of positron scattering and annihilation on noble-gas 
atoms. In the future it will be important to calculate explicitly the contribution of 
the RPA-type screening diagrams and also to move to a fully relativistic framework, 
particularly for heavy atoms and ions. 
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